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Abstract We present the outcome of a highly-sensitive search for magnetic fields 
on the cool supergiant Betelgeuse. A time-series of six circularly-polarized spectra 
was obtained using the NARVAL spectropolarimeter at Telescope Bernard Lyot (Pic 
du Midi Observatory), between 2010 March and April. Zeeman signatures were re- 
peatedly detected in cross-correlation profiles, corresponding to a longitudinal com- 
ponent of about 1 G. The time-series unveils a smooth increase of the longitudinal 
field from 0.5 to 1.5 G, correlated with radial velocity fluctuations. We observe a 
strong asymmetry of Stokes V signatures, also varying in correlation with the radial 
velocity. The Stokes V line profiles are red-shifted by about 9 kms^' with respect 
to the Stokes 1 profiles, suggesting that the observed magnetic elements may be 
concentrated in the sinking components of the convective flows. 
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Fig. 1 Time-series of Stokes V (bold lines) and their corresponding "null" line profiles (thin line 
below). The successive profiles are shifted vertically for display clarity. After Auriere et al. (2010). 

1 Introduction 

The widespread signatures of magnetic activity in cool stars are believed to origi- 
nate from the coexistence of convection and rotation in their outer layers. Since the 
first dynamo models of Parker (1955), rotation is generally accepted as a central 
ingredient of these stellar dynamos, in a two-level action. The stellar spin is first in- 
volved in the generation of radial and latitudinal shears that are able to wind up the 
field lines of a seed magnetic field around the rotation axis, resulting in the creation 
of a strong toroidal field component. Stellar rotation is also acting on the vertical 
plasma flows through the Coriolis force that succeeds at generating helical motions 
able to twist again the field lines of the toroidal field and generate a poloidal field 
component. 

If many details of the physical processes involved in such large-scale dynamos 
are still a matter of debate (see e.g. Charbonneau 2010 for a review), this theoreti- 
cal framework is now widely accepted to interpret the cycUcal activity behaviour of 
the Sun and solar-type stars. However, other models suggest that convection alone 
is able to sustain a dynamo, without any rotational effects involved (e.g. Cattaneo 
1999, Vogler & Schiissler 2007). Turbulent dynamo action may be responsible for 
the smallest-scale magnetic elements observed on the solar surface (Lites et al. 
2008), although intranetwork magnetic elements may also result from the decay 
of larger magnetic regions (created by the global dynamo) through the continuous 
convective mixing of the solar upper layers. 
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Since both global and small-scale dynamos may be simultaneously active in the 
Sun, it is not easy to disentangle the respective magnetic outcome of these two 
different processes. A promising way to reach this goal consists in observing a star 
with no rotation at all, or at least a star rotating so slowly that the onset of a global 
dynamo in its internal layers is unlikely. If stellar spectropolarimetry is our best 
asset to detect a magnetic field in a non-rotating star, the polarimetric detection of 
Zeeman signatures is mostly insensitive to small-scale magnetic elements as those 
expected to be generated by a local dynamo. This issue is unescapable for solar-type 
dwarfs, on which millions of photospheric convective cells are visible at any time, 
resulting in a highly tangled intranetwork field pattern. Cool supergiant stars may 
offer a rare opportunity to circumvent this problem, since their convective cells are 
expected to be much larger than on the Sun, with only a few of them covering the 
stellar surface (Schwarzschild 1975, Chiavassa et al. 2010), so that the spatial scale 
of convection may be sufficiently large to limit the mutual cancellation of Zeeman 
signatures of close-by magnetic elements with opposite polarities. 

To test the feasabiUty of magnetic field detection in cool supergiant stars, we 
have concentrated on Betelgeuse (a Orionis), one of the brightest members of this 
class. We briefly summarize a few well-known properties of this object, present 
the spectropolarimetric observations gathered with NARVAL and discuss the first 
outcome of this project. 



2 The cool supergiant a Orionis 

With a M2Iab spectral type classification, Betelgeuse can be taken as the prototype 
of cool supergiant stars. Thanks to its proximity to the Earth, it was the first star 
to have its radius determined using interferometry (Michelson & Pease 1921), and 
a recent value of 645 ± 129/?0 was derived in near-infrared interferometry (Perrin 
et al. 2004). A significant scatter is observed in radius estimates depending on the 
adopted wavelength domain, because of the extension and thermal structure of the 
stellar atmosphere (see e.g. Uitenbroek et al. 1998, and references therein). The 
efficient mass-loss of Betelgeuse is at the origin of an inhomogeneous distribution 
of gas and dust, detected up to a few tens of stellar radii (Kervella et al. 201 1). 

Using spatially-resolved, high-resolution UV spectroscopy with the HST, Uiten- 
broek et al. (1998) were able to propose a rotation period of about 17 yr and a low 
incUnation of the rotation axis, of about 20°. The basic ingredients of a global dy- 
namo, primarily based on rotation, are therefore likely missing in a Orionis. 
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Fig. 2 Estimates of the disc-integrated, line-of-sight component of the magnetic field, from 2010 
March 14 to 21010 April 17 (Auriere et al. 2010). 



3 Surface magnetic field 

Spectropolarimetric observations of Betelgeuse were obtained using the NARVAL 
spectropolarimeter (Auriere 2003), installed at Telescope Bernard Lyoj^ Pic du 
Midi Observatory (France), in a highly-sensitive search for a weak surface mag- 
netic field (Auriere et al. 2010). The data were collected during 6 nights in 2010, 
from March 14 to April 17. They consist of high-resolution (R ~ 65,000) echelle 
spectra in light intensity (Stokes I) and circular polarization (Stokes V), offering 
an almost contiguous coverage of the wavelength interval between 370 nm and 
1,000 nm. Each Stokes V spectrum is built from a sequence of 4 exposures taken 
with different azimuths of the polarimetric optics. Using a different combination of 
the sub-exposures, a "null" control spectrum is also computed, which should not 
contain any detectable signature. To reach a high signal-to-noise ratio and avoid any 
saturation of the detector, a sequence of 16 to 20 spectra was acquired for every 
individual night and later averaged. 

To lower further the noise level and improve the detectability of Zeeman sig- 
natures, all spectra were processed using the Least-Squares-Deconvolution (LSD) 
cross-correlation technique (Donati et al. 1997, Kochukhov et al. 2010). By doing 
so, mean line profiles were computed, using a list of 15,000 photospheric atomic 
lines corresponding to the atmospheric parameters of Betelgeuse. The resulting 
time-series of Stokes V profiles is plotted in Fig.[T| Circularly polarized signatures. 



The Bernard Lyot Telescope is operated by the Institut National des Sciences de I'Univers of the 
Centre National de la Recherche Scientifique of France. 
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most likely generated through the Zeeman effect, are detected above noise level at 
the radial velocity of the star. 

Using the centre-of-gravity technique (Rees & Semel 1979), the Stokes V sig- 
natures can be translated into estimates of the line-of-sight component (B/) of the 
magnetic field (Auriere et al 2010). The series of B/ measurements are plotted in 
Fig. |2] They display an average value of the order of 1 G. If the variability of Stokes 
V profiles is barely visible to the naked eye, Bi estimates enable one to unveil a 
first type of variability of the Zeeman signatures, with a regular increase of the field 
strength over our observing window, from 0.5 to 1.5 G. 



4 Surface convection and magnetic activity 

The wide spectral coverage and high spectral resolution of NARVAL spectra give 
access to a number of classical activity indicators that provide us with a useful set 
of additional measurements, carrying information that complement the longitudinal 
field estimates. 



4.1 Chromospheric emission 

The stellar chromospheric flux can be inferred from the various chromospheric lines 
showing up in NARVAL intensity spectra. We choose here to concentrate on the Ca 
II infrared triplet, located in a spectral domain where the signal-to-noise ratio is 
generally high in Betelgeuse spectra. 

An emission index is constructed by integrating the fluxes Ca\, Ca2 and Cut, 
in three rectangular bandpasses, centered around the three component of the Ca II 
triplet (at 849.8, 854.2 and 866.2 nm), each one with a width of 0.2 nm. We also 
integrate the flux in two rectangular bands Ci and C2 with a width of 0.5 nm in the 
neighbouring spectral domain (centered around 847.58 and 870.49 nm) to normalize 
the index, which is therefore computed as follows: 

C1+C2 

The temporal evolution of this index is plotted in Fig. [3] Significant fluctuations are 
observed. If a global decrease is visible between March and April, the evolution is 
not as smooth as the increase of B/. The Ca II emission is generally expected to be 
correlated to the local magnetic field strength, so that the (loose) anti-correlation 
between B/ and the Ca II index might seem paradoxical. We stress, however, that in 
the case of stars with low vsin/ values, the polarimetric signatures are specifically 
sensitive to the low-order component of the magnetic topology, while the chromo- 
spheric emission retains information about a wide range of spatial scales (Petit et 
al. 2008). The different temporal trend of both tracers is therefore not a surprise in 
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Fig. 3 Additional activity indicators derived from NARVAL spectra. From top to bottom, we plot 
the Ca II IRT index, the mean radial velocity of Stokes I LSD profiles (RV) in km s^ ' , the velocity 
span of bisectors calculated from Stokes I LSD profiles (kms^'), the relative surface asymmetry 
(dA) and the relative amplitude asymmetry (da) of Stokes V LSD profiles. Statistical error bars are 
not indicated in the plots whenever they are smaller than the symbol size. 



the context of a magnetic field probably shaped by a complex pattern of convective 
flows. 



4.2 Radial velocity and profile bisectors 

The radial velocities (RV) of Stokes I LSD profiles, derived from the fitting of a 
gaussian on the line profiles, are plotted in Fig. [3] following previous work by Gray 
(2008). The wavelength calibration is performed using telluric lines recorded in 
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NARVAL spectra. For solar-type dwarfs, the resulting RV stability is of the order 
of 15 to 30 ms^' (Moutou et al. 2007). A similar accuracy can be expected for 
Betelgeuse, for which LSD profiles benefit from a much larger number of spectral 
lines. The observed fluctuations are far above this limit, reaching about 1.5 kms^'. 
The evolution is smooth and anti -correlated to B/ measurements. 

The Stokes I LSD profiles of Betelgeuse are also highly asymmetric, with a red 
wing of the line much wider than the blue wing (Auriere et al. 2010). Bisectors 
computed from the LSD profiles confirm this trend, with a clear deviation towards 
the red near the continuum. The total velocity span of the bisectors (plotted in Fig. 
[sjl is comprised between 1 and 2 kms^^ Fluctuations are less organized than those 
recorded in RV, with a marginal trend to display a lower bisector span in April. 



4.3 Stokes V asymmetry and zero-crossing of Stokes V profiles 

The double-peaked Stokes V profiles of Betelgeuse all possess a significant level of 
asymmetry, in the sense that the blue and red lobes of the profiles have a different 
amplitude and delimit a different surface. This property of Stokes V signatures is 
widely observed on the Sun (Viticchie & Sanchez Almeida 2011) and is generally 
linked to vertical gradients of both velocity and magnetic field in magnetic elements 
(Lopez Ariste 2002). 

To quantify this observation, we use the approach of Petit et al. (2005) and call 
Uh and Or the amplitude of the blue and red lobes of the profile. We then derive 
a relative amplitude asymmetry da — [at, — ar)/{ab + Or). In a similar manner, we 
define a relative surface asymmetry dA = {At, —Ar)/{Ah +Ar), where A/, and Ar are 
the areas of the blue and red lobe, respectively. The successive values of dA and da 
are plotted in Fig. [3] We observe a global decrease of both parameters, with a more 
regular trend in da. Both parameters change sign during the time-series. 

The double-lobed shape of Stokes V profiles implies a sign reversal close to the 
line-core. The zero-crossing radial velocity observed for the Stokes V profiles of 
Betelgeuse is of about 33 kms^' (with no clear temporal trend), while the core of 
the Stokes I profile is located at a significantly smaller velocity of about 24 kms^'. 
These systematic velocity shifts between Stokes I and V can have different origins. 
A first possible explanation would involve the Stokes V asymmetry reported above, 
as it can induce a shift of the zero-crossing in the case of a limited spectral resolution 
(Solanki & Stenflo 1986); but in this case, the sharp evolution of da and dA should 
progressively displace the zero-crossing towards lower radial velocities, which is 
not observed. A second option involves a large-scale toroidal magnetic field (Petit 
et al. 2005), but the existence of a toroidal component would be hard to reconcile 
with the absence of any significant rotation in Betelgeuse (unless the slowly rotat- 
ing surface is hiding faster rotation in internal layers). As a third option, we propose 
that the red-shift of Stokes V profiles may be related to the concentration of mag- 
netic elements in the sinking component of the convective mixing, as observed in 
small-scale solar magnetic elements. If this last interpretation seems more easy to 
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reconcile with the basic properties of Betelgeuse and is also quantitatively consis- 
tent with the convective velocity measurements of Gray (2008), further investigation 
is obviously needed before reaching any conclusion. 



5 A small-scale dynamo in Betelgeuse ? 

In the past, strong magnetic fields have been detected in fast-rotating giants or sub- 
giants belonging to RS CVn systems or to the FK Com class (Petit et al. 2004a, 
2004b). More recently, repeated magnetic field detections have been obtained in 
active, single red giants (Konstantinova et al. 2008, 2010, Auriere et al. 2009, 
Sennhauser & Berdyugina 2011) and in a Ukely descendant of a strongly magnetic 
Ap star (Auriere et al. 2008). 

In this growing literature on evolved stars, the detection of a weak surface mag- 
netic field at the surface of Betelgeuse (Auriere et al. 2010) is an important observa- 
tional result, in the sense that the physical interpretations proposed for other objects 
to account for their magnetic nature cannot be applied here. Firstly, the magnetic 
field of Betelgeuse has to be generated without the help of a fast, or even moderate 
stellar rotation, and this specificity should exclude any global dynamo. Secondly, 
the very large radius implies that any magnetic remnant of a strong magnetic field 
on the main sequence would be too diluted to be detectable at photospheric level. In 
this situation, a more natural interpretation would involve the convection alone as 
the engine of a dynamo, bringing the first strong observational evidence that such 
process (proposed by Dorch & Freytag 2003) can be efficient in cool stars. 

This exciting result, later confirmed for a larger sample of cool supergiants 
(Grunhut et al. 2010), comes together with a number of additional tracers of mag- 
netic activity and convection (chromospheric emission, radial velocities, line bisec- 
tors. Stokes V asymmetries). This wealth of information is a motivation to pursue 
the spectropolarimetric monitoring of Betelgeuse, in order to investigate longer-term 
trends that may affect the various measurements at our disposal and study the pos- 
sible role of the surface magnetic field in the onset of the mass-loss of Betelgeuse 
and other supergiant stars. 



References 

1. Auriere M., 2003, EAS, 9, 105 

2. Auriere M., et al., 2008, A&A, 491, 499 

3. Auriere M., et al., 2009, A&A, 504, 231 

4. Auriere M., Donati J.-R, Konstantinova-Antova R., Perrin G., Petit P., Roudier T., 2010, 

A&A, 516, L2 

5. Cattaneo R, 1999, ApJ, 5 15, L39 

6. Charbonneau R, 2010, LRSP, 7, 3 

7. Chiavassa A., Haubois X., Young J. S., Plez B., Josselin E., Perrin G., Freytag B., 2010, A&A, 
515, A12 



Magnetic field and convection in the cool supergiant Betelgeuse 



9 



8. Donati J.-R, Semel M., Carter B. D., Rees D. E., Collier Cameron A., 1997, MNRAS, 291, 
658 

9. Dorch S. B. R, Freytag B., 2003, lAUS, 210, 12P 

10. Gray D. R, 2008, AJ, 135, 1450 

11. Grunhut J. H., Wade G. A., Hanes D. A., Alecian E., 2010, MNRAS, 408, 2290 

12. Kervella P., Perrin G., Chiavassa A., Ridgway S. T., Cami J., Haubois X., Verhoelst T., 2011, 

A&A, 531, A117 

13. Kochukhov O., Makaganiuk V., Piskunov N., 2010, A&A, 524, A5 

14. Konstantinova-Antova R., Auriere M., lUev I. K., Cabanac R., Donati J.-F., Mouillet D., Petit 
P., 2008, A&A, 480, 475 

15. Konstantinova-Antova R., et al., 2010, A&A, 524, A57 

16. Lites B. W., et al., 2008, ApJ, 672, 1237 

17. Lopez Ariste A., 2002, ApJ, 564, 379 

18. Michelson A. A., Pease R G., 1921, ApJ, 53, 249 

19. Moutou C, et al., 2007, A&A, 473, 651 

20. Parker E. N., 1955, ApJ, 122, 293 

21. Perrin G., Ridgway S. T., Coude du Foresto V., Mennesson B., Traub W. A., Lacasse M. G., 
2004, A&A, 418, 675 

22. Petit R, et al., 2004, MNRAS, 348, 1175 

23. Petit R, et al, 2004, MNRAS, 351, 826 

24. Petit P, et al., 2005, MNRAS, 361, 837 

25. Petit P, et al., 2008, MNRAS, 388, 80 

26. Rees D. E., Semel M. D., 1979, A&A, 74, 1 

27. Schwarzschild M., 1975, ApJ, 195, 137 

28. Sennhauser C, Berdyugina S. V., 2011, A&A, 529, AlOO 

29. Solanki S. K., Stenflo J. O., 1986, A&A, 170, 120 

30. Uitenbroek H., Dupree A. K., Gilliland R. L., 1998, AJ, 116, 2501 

31. Viticchi6 B., Sdnchez Almeida J., 2011, A&A, 530, A14 

32. Vogler A., Schussler M., 2007, A&A, 465, L43 



